Coupled Chemistry- Emission Model for Atomic Oxygen Green 
and Red-doublet Emissions in Comet C/1996 B2 Hyakutake 
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ABSTRACT 

The green (5577 A) and red-doublet (6300, 6364 A) lines are prompt emissions of metastable 
oxygen atoms in the *S and 1 D states, respectively, that have been observed in several comets. 
The value of intensity ratio of green to red-doublet (G/R ratio) of 0.1 has been used as a bench- 
mark to identify the parent molecule of oxygen lines as H 2 0. A coupled chemistry-emission model 
is developed to study the production and loss mechanisms of O^S) and 0( X D) atoms and the 
generation of red and green lines in the coma of C/1996 B2 Hyakutake. The G/R ratio depends 
not only on photochemistry, but also on the projected area observed for cometary coma, which 
is a function of the dimension of the slit used and geocentric distance of the comet. Calculations 
show that the contribution of photodissociation of H 2 to the green (red) line emission is 30 
to 70% (60 to 90%), while CO2 and CO are the next potential sources contributing 25 to 50% 
(<5%). The ratio of the photo-production rate of 0( 1 S) to 0( 1 D) would be around 0.03 (± 
0.01) if H2O is the main source of oxygen lines, whereas it is ^0.6 if the parent is CO2. Our 
calculations suggest that the yield of O^S) production in the photodissociation of H2O cannot 
be larger than 1%. The model calculated radial brightness profiles of the red and green lines and 
G/R ratios are in good agreement with the observations made on comet Hyakutake in March 
1996. 

Subject headings: comets: general - comets: individual (C/1996 B2 Hyakutake) - Molecular processes 



1. Introduction 

The spectroscopic emissions from dissociative 
products in cometary coma are often used in es- 
timating production rates of respective cometary 
parent species which are sublimating directly from 



the nucleus (Feldman et al. 2004 Combi et al. 
2004). It is a known fact that at smaller (<2 



it's abundance indirectly based on the emissions 
from daughter products, like OH, O and H. Thus, 
tracking emissions of the dissociative products of 
H2O has became an important diagnostic tool in 
estimating the production rate as well as in under- 
standing the spatial distribution of H 2 in comets 



AU) heliocentric distances, the inner cometary 
coma is dominantly composed of H 2 0. The in- 
frared emissions of H2O molecule are inaccessible 
from ground because of strong attenuation by the 
terrestrial atmosphere. Since H2O does not show 
any spectroscopic transitions in ultraviolet or vis- 
ible regions of solar spectrum, one can estimate 



pelsemme fc Combi||1976[|1979[|Fink fc Johnson 



1984 Sc hultz et al.|1992||Morgenthaler et al.|2001| 



Furusho et al. 2006). For estimating the density 
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distribution of H2O from the emissions of daugh- 
ter species, one has to account for photochemistry 
and associated emission processes. 

The major dissociative channel of H2O is the 
formation of H and OH, but a small fraction is also 
possible in 0( 3 P, 1 S, 1 D) and H 2 . The radiative 
decay of metastable 1 D and *S states of atomic 



oxygen leads to emissions at wavelengths 6300, 
6364 A (red doublet) and 5577 A (green line), 
respectively. The energy levels of atomic oxygen 
and these forbidden transitions are shown in Fig- 
ure [l] Even though these emissions are accessi- 
ble from ground-based observatories, most of the 
times they are contaminated by telluric night sky 
emissions as well as emissions from other cometary 
species. Doppler shift of these lines, which is a 
function of the relative velocity of comet with re- 
spect to the Earth, offers a separation from telluric 
emissions provided a high resolution cometary 
spectrum is obtained. In most of the cometary ob- 
servations it is very difficult to separate the green 
line in optical spectrum because of the contamina- 
tion from cometary C2 (1-2) P-branch band emis- 
sion. The red line 6300 A emission is also mildly 
contaminated by the Q-branch emission of NH2 
molecule, but in high resolution spectrum this can 
be easily resolved. 

Since these atomic oxygen emissions result due 
to electronic transitions which are forbidden by se- 
lection rules, solar radiation cannot populate these 
excited states directly from the ground state via 
resonance fluorescence. The photodissociative ex- 
citation and electron impact excitation of neu- 
tral species containing atomic oxygen, and ion- 
electron dissociative recombination of O-bearing 
i on species, can produce the se metastable states 
([Bhardwaj fc Haider| [2002]) . If 0( 1 D) is not 



quenched by ambient cometary species, then pho- 
tons at wavelengths 6300 and 6364 A will be emit- 
ted in radiative decay to the ground 3 P state. 
Only about 5% of 0( 1 S) atoms result in 2972 and 
2958 A emissions via direct radiative transition to 
the ground 3 P state of atomic oxygen. Around 
95% of 0( 1 S) decays to the ground state through 
O^D) by emitting green line (cf. Fig. nil. This 
implies that if the green line emission is present 
in cometary coma, the red doublet emission will 
also be present, but the opposite is not always 
true. The average lifetime of 0( X D) is relatively 
small (~110 s) compared to the lifetime of H2O 
molecule (-8 x 10 4 s) at 1 AU. The 0( X S) also has 
a very short average lifetime of about 0.1 s. Due to 
the short lifetime of these metastable species, they 
cannot travel larger distances in cometary coma 
before de-exciting via radiative transitions. Hence, 
these emissions have been used as diagnostic tools 



to estimate the abundance of H2O in comets (Fink 



fc Johnson|1984|IMagee-Sauer et al.|1990||Morgen- 
thaler et al.|2001 ). The intensity of O [I] emissions, 
in Rayleigh, can be calculated using the following 
equation ( Festou fc Feldmaii||1981 ) 



/ = lO^T^afiN 



(1) 



where t p is the lifetime of excited species in sec- 
onds, a is the yield of photodissociation, /3 is the 
branching ratio, and N is the column density of 
cometary species in cm~ 2 . 

In the case of red doublet (6300 and 6364 A), 
since both emissions arise due to transition from 
the same excited state (2P 4 1 D) to the ground 
triplet state (2P 4 3 P), the intensity ratio of these 
two lines should be the same as that of branching 
ratio of corresponding transitions. Using Einstein 



transition probabilities, Storey & Zeippen (2000) 



calculated the intensity ratio of red doublet and 
suggested that the intensity of 6300 A emission 
would be 3 times stronger than that of 6364 A 
emission, and this has been observed in several 



comets also ( Spinradl 1982 Fink & Johnson 1984 



Mor rison et al.||1997| |Cochran fc Cochran||2001 



Capria et al.j 2005; Furusho et al.||2006 
et al.||2008t |Cochran|2008] ) 



Capria 



The ratio of intensity of green line to the sum 
of intensities of red doublet can be calculated as 



'5577 



I, 



6300 



'6364 



T green Qreen^* greenP green 
T red a redN re d(f36300+6364) 



(2) 



If the emission intensities of oxygen lines are com- 
pletely attributed to only photodissociative excita- 
tion of H2O and column densities are assumed al- 
most same for both emissions, then the ratio of in- 
tensities of green line to red dou blet is directly pro- 
portional to the ratio of T~ l aj3. Festou &; Feldman 



( 1981 ) reviewed these atomic oxygen emissions in 
comets. Based on the observation of 0[I] 2972 A 
emission in the IUE spectrograph of comet Brad- 



field (1979X), [Festou fc Feldnrah] ( |1981[ ) calculated 
the brightness profiles of red and green emissions. 



Festou & Feldman (1981) also calculated a theo- 



retical value for the ratio of the intensity of green 
line to red doublet (hereafter refer to as the G/R 
ratio), which has a value of around 0.1 if H2O is 
the source for these 0[I] emissions in cometary co- 
mae, and it is nearly 1 if the source is CO2 or 
CO. Observations of green and red line emissions 
in several comets have shown that the G/R ratio 



is around 0.1, suggesting that H 2 is the main 
source of these 0[I] lines. However, since no ex- 
perimental cross section or yield for the produc- 
tion of Of^S) from H2O is available in literature, 



the G/R ratio has been questioned by Huestis & 
Slanged ([2006]). 



Generally, the red line is more intense than the 
green line because the production of 0( 1 D) via 
dissociative excitation of H2O is larger compared 
to the radiative decay of 0( 1 S). Since the lifetime 
of 0( 1 D) is larger, quenching is also a significant 
loss process near the nucleus. So far, the observed 
G/R ratio in comets is found to vary from 0.022 



to 0.3 (Cochran 1984, 2008 Morrison et al. 1997 



Zhang et al ||2001[ |Cochran fc Cochran||200i| |Fu 



I 



rusho et al.|2006||Capria et al.||2005[ |2008[ |2010| ) 

There are several reactions not involving H2O 
which can also produce these forbidden oxygen 
lines (Bhardwaj & Haider [2002). Among the O- 



bearing species, CO2 and CO also have dissocia- 
tive channels producing O^D) and 0( 1 S). How- 
ever, complex O-bearing molecules (e.g., H 2 CO, 
CH3OH, HCOOH) do not produce atomic oxy- 
gen as a first dissociative product. Based on the 



brightness of 6300 A emission intensity, Delsemme 



& Combi 



(1976) derived the production rate of 
0( i D) in comet Bennett 1970 II and suggested 
that the abundance of CO2 is more than that of 



H 2 0. Delsemme & Combi (1979) estimated the 
production of 0( i D) in dissociation of H2O and 
CO2; about 12% of H2O is dissociated into H2 and 
0( 1 D), while 67% of C0 2 is dissociated into CO 
and O^D). They suggested that a small amount 
of CO2 can contribute much more than H 2 to 
the red doublet emission. The model calculations 



of Bhardwaj & Haider ( 2002 ) showed that the pro- 
duction of 0( 1 D) is largely through photodissocia- 
tive excitation of H2O while the major loss mecha- 
nism in the innermost coma is quenching by H2O. 



Cochran & Cochran (2001 ), based on the observa- 



tion of width of red and green lines, argued that 
there must be another potential source of atomic 
oxygen in addition to H2O, which can produce 
0( 1 S) and O^D). Observations of the green and 
red lines in nine comets showed that the green line 
is wider than the red line (Cochran 20081, which 



could be because various parent sources are in- 
volved in the production of 0( 1 S). 

The model of |Glinski et al.| ( |2004[ ) showed 
that the chemistry in the inner coma can pro- 



duce 1% O2, which can also be a source of red 
and green lines. Manfroid et al. (20071 also ar- 



gued, based on lightcurves, that forbidden 0[I] 
emissions are probably contributed through dis- 
sociation sequence of CO2. Recent observation of 
comet 17P/Holmes showed that the G/R ratio can 
be even 0.3, which is the highest reported value 
so far: suggesting that C0 2 and CO abundances 



might be higher at the time of observation ( Capria 
eTaLl[2010| . 



Considering various arguments based on differ- 
ent observations and theoretical works, we have 
developed a coupled chemistry-emission model 
to quantify various mechanisms involved in the 
production of red and green line emissions of 
atomic oxygen. We have calculated the produc- 
tion and loss rates, and the density profiles, of 
metastable O(D) and O^S) atoms from the O- 
bearing species, like H 2 0, C0 2 , and CO, and also 
from the dissociated products OH and O. This 
model is applied to comet C/1996 B2 Hyakutake, 
which was studied through several observations in 



1999 



1996 March (|Biver et al 

T9971 ICochran & Cochranl [2001 



Morrison et al. 



Morgenthaler 



et al.|2001[|Combi et al.l2005|[Cochran|2008[ ). The 
line-of-sight integrated brightness profiles along 
cometocentric distances are calculated for 5577 
and 6300 A emissions and compared with the ob- 
served profiles of Cochran ( 2008| ). We have also 
evaluated the role of slit dimension, used in the 
observation, in determining the G/R ratio. The 
aim of this study is to understand the processes 
that determine the value of G/R ratio. 

2. Model 

The neutral parent species considered in this 
model are H 2 0, CO2, and CO. We do not consider 
other significant O-bearing species, like H2CO, 
CH3OH, since their first dissociation does not lead 
to the formation of atomic oxygen atom; the O 
atom appears in subsequent photodissociation of 
daughter products, like OH, CO, HCO. On 1996 
March 24, the H2O production rate for comet 



C/199 6 B2 Hyakutake measured by Mumma et al. 



(1996) was 1.7 x 10^ y s 



Based on H Ly-a 



emission observation, Combi et al 



sured H2O production rate as 2.6 



(1998) 
10^s 



mea- 
-1 ™ 



1996 April 4. Using molecular radio line emissions, 



Biver et al. ( 1999 ) derived the production rates of 



different species at various heliocentric distances 
from 1.6 to 0.3 AU. They found that around 1 AU 
the relative abundance of CO with respect to H 2 
is high (-22%) in the comet C/1996 B2 Hyaku- 
take. 

The number density n^(r) of i parent species 
at a cometocentric distance r in the coma is cal- 
culated using the following Hascr's formula 



rii(r) 



47ru,-r' 



■(e'^ r ) 



(3) 



Here Q p is the total gas production rate of the 
comet, Vi and fy are the gas expansion velocity 
(taken as 0.8 

5.0 



km s 
scale length (/3h 2 o 



Biver et al.| 1999) and the 
10 4 km, /3co 2 



8.2 x 
= 1.4 



10 6 km) of the 



lh 



x 10 5 km, and (3 C o 
species, respectively. The Haser model's neutral 
density distribution has been used in several pre- 
vious studies for deriving the production rate of 
H 2 in comets based on the intensity of 6300 A 



emission (Delsemme & Combi 1976 



& Johnson 1984 Morgenthaler et al. 2001). In 



1979 Fink 



our model calculations the H2O production rate 
on 1996 March 30 is taken as 2.2 x 10 29 s" 1 . The 
abundance of CO relative to H 2 is taken as 22%. 
Since there is no report on the observation of C0 2 
in the comet Hyakutake, we assumed its abun- 
dance as 1% relative to H 2 0. However, we vary 
C0 2 abundance to evaluate its effect on the green 
and red-doublet emissions. The calculations are 
made when the comet C/1996 B2 Hyakutake was 
at a heliocentric distance of 0.94 AU and a geocen- 
tric distance of 0.19 AU on 1996 March 30. The 
calculated G/R ratio on other days of the obser- 
vation is also reported. 

The number density of OH produced in disso- 
ciation of parent species H 2 at a given cometo- 
centric distance r is calculated using Haser's two 
parameter coma model 



1 s Qp Pp , _ 

n H(r) = 2 « ^( e 



ppr 



-Pnr\ 



(4) 



Here v is the average velocity of daughter species 
taken as 1 km s" 1 , and /3p and /3r are the destruc- 
tion scale lengths of the parent (H 2 0, 8.2 x 10 4 
km) and daughter (OH, 1.32 x 10 5 km) species, 
respectively ( Huebner et al.|[l992 ). The solar UV- 
EUV flux is taken from SOLAR 2000 v.2.3.6 (S2K) 
model of iTobiska et all (12000} for the day 1996 



March 30, which is shown i n Figure [2} Fo r com- 
parison the solar flux used by Huebner et al. ( 1992 ) 
in calculating 0( 1 D) and 0( i S) production rates 
from various O-bearing species is also presented in 
the same Figure. 

The primary photoelectron energy spectrum 
Q(E,r,9) is calculated by degrading solar radia- 
tion in the neutral atmosphere using 

Q(E, r,0)=J2 [ "*M ^ A ) /o °( A ) ex Pi- T ( r > e > A )l dX 

(5) 
where, 

ni {r')dr' (6) 

Here crf(X) and erf (A) are the absorption and 
ionization cross sections, respectively, of the i th 
species at the wavelength A, nj(r) is its neutral 
gas density and r(r, 9, A) is optical depth of the 
medium at the solar zenith angle 6. /00(A) is the 
unattenuated solar flux at the top of atmosphere 
at wavelength A. All calculations are made at so- 
lar zenith angle 9 of 0°. The total photoabsorption 
and photoionization cross sections of H 2 0, C0 2 , 
and CO are taken from the compilation of Hueb- 



ner et al. ( 1992 ) (http : //amop . space . swri . edu) 



and interpolated at 10 A bins to make them com- 
patible with the S2K solar flux wavelength bins 
for use in our model calculations. The total pho- 
toabsorption and photoionization cross sections 
for H 2 0, C0 2 , and CO are presented in Figure [3] 
The photochemical production rates for ionization 
and excitation of various species are calculated us- 
ing degraded solar flux and cross sections of cor- 
responding processes (discussed in Section [3]) at 
different cometocentric distances. 

The primary photoelectrons are degraded in 
cometary coma to calculate the steady state pho- 
toelectron flux using the Analytical Yield Spec- 
trum (AYS) approach, which is based on the 
Monte Carlo method ([Singhal k Bhardwaj||1991 
Bhardwaj & Singhal 1993; Bhardwaj & Michae 



1999a| |Bhardwaj fc Jain||2009[ ) . The AYS method 
of degrading electrons in the neutral atmosphere 
can be explained briefly in the following manner. 
Monoenergetic electrons incident along Z-axis in 
an infinite medium are degraded in collision-by- 
collision manner using the Monte Carlo technique. 
The energy and position of the primary electron 



and its secondary or tertiary are recorded at the 
instant of an inelastic collision. The total num- 
ber of inelastic events in the spatial and energy 
bins, after the incident electron and all its sec- 
ondaries and tertiaries have been completely de- 
graded, is used to generate numerical yield spec- 
tra. These yield spectra contain the yield informa- 
tion about the electron degradation process and 
can be employed to calculate the yield for any in- 
elastic event. The numerical yield spectra gener- 
ated in this way are in turn represented analyt- 
ically, which contains the information about all 
possible collisional events based on the input elec- 
tron impact cross sections, resulting in the AYS. 
This yield spectrum can be used to calculate the 
steady state photoelectron flux. More details of 
the AYS approach and the method of photoelec- 
tron computation are given in several previous 



1990 



1999 



fc B 



papers (Singhal & Haider 1984| Bhardwaj et al 



I996| JSinghal fe Bhardwaj||1991| |Bhardwaj 
2003] |Bhardwaj fc Michael||1999bi IHaiderl 



hardwaj 2005; Bhardwaj fc Jain 2009 2011 



Raghuram fc Bhardwaj|2011 ) . The total inelastic 



electron impact cross sections for H 2 are taken 
from |Jackman et al. (19771 and Seng & Linder 
( 1976 ), and those for CO2 and CO are taken from 



Jackman et al. (1977). The electron impact cross 



sections for different dissociative ionization states 



of H2O are taken from Itikawa & Mason (20051, 



for C0 2 from Bhardwaj & Jain (2009), and for 



CO from McConkey et al. (2008). The volume 



excitation rates for different processes are calcu- 
lated using steady state photoelectron flux and 
electron impact cross sections. The electron tem- 
perature required for ion-electron dissociative re- 
combination reactions is taken from |Korosmezey| 
et al. (1987). The detailed description of cou- 



pled chemistry-transport model has been given in 



our earlier papers (Bhardwaj et al. 1995 



Bhardwaj||1999l|Bhardwaj fc Haider||2002 



fc Bhardwaj ||2005{ |Bhardwaj fc Raghuram 



1996 



Haider| 



2011) 



Various reactions involved in the production and 
loss of metastable 0( 1 S) and 0( 1 D) atoms con- 
sidered in our model are listed in Tables [l] and [2] 
respectively. 



3. Dissociation of neutral species produc- 
ing 0( 1 S) and O^D) 

3.1. Photodissociation 

3.1.1. H 2 and OH 

The dissociation of H 2 molecule starts at 
wavelengths less than 2424 A and the primary 
products are H and OH. But the pr e-dissociation 
process mainly starts from 1860 A (Watanabe & 



Zelikoff 1953). The threshold wavelength for the 
photoionization of H 2 is 984 A. Hence, solar UV 
photons in the wavelength region 1860 to 984 A 
can dissociate H 2 and produce different daugh- 
ter products. The threshold wavelengths for the 
dissociation of H 2 resulting in the production of 
0( 1 S) and 0( 1 D) are 1390 A and 1770 A, respec- 
tively. Till now, the photo-yield value for the pro- 
duction of 0( 1 D) from H 2 have been measured 



in only two experiments. Slanger & Black ( 1982 ) 



measured the 0( 1 D) yield in photodissociation of 
H 2 at 1216 A, and found its value to be 10%. 



McNesby et al.| ([1962) reported a 25% yield for 
the production of OfT)) or 0( 1 S) at 1236 A from 
H 2 0. 



Huebner et al. (19921 calculated photo pro- 



duction rates for different excited species pro- 
duced from H 2 using absorption and ionization 
cross sections compiled from different experimen- 
tal measurements. In our model the cross sec- 
tions for the production of 0( 1 D) in photodis- 
sociation of H 2 are taken from |Huebner et al. 
(1992), which were determined based on experi- 
ments of |Slanger fc Black| (|1982[) and |McNesby 
et al. (1962). Huebner et al. ( 1992[ ) assumed that 
in the 1770 to 1300 A wavelength region around 
25% of H 2 molecules photodissociate into H 2 and 
0( 1 D), while between 1300 and 984 A about 10% 
of H 2 Q disso ciation produces Of 1 ! )) (cf. Fig. [4J>. 
Below 984 A, Huebner et al. ( 1992 ) assumed that 



33% of di ssociation of H 2 Q leads to the formation 
of 0( 1 D). Festou (1981) discussed various dissoci- 



ation channels for H 2 in the wavelength region 
less than 1860 A. Solar photons in the wavelength 
region 1357 to 1860 A dissociates around 72% of 
H 2 molecules into ground states of H and OH. 



But, according to Stief et al. ( 1975 ) approximately 



1% of H 2 molecules are dissociated into H 2 and 
0( 1 D) in this wavelength region. The calculated 
rates for the 0( 1 D) production from photodisso- 



dative excitation of H 2 by Huebner et al. ( 1992 ) 



are 5.97 x 1CT 7 s^ 1 and 1.48 X 10 _e s" 1 for so- 
lar quiet and active conditions, respectively. Using 
the S2K solar EUV-UV flux on 1996 March 30 and 



cross sections from Huebner et al. ( 1992 ) (see Fig- 



ure |4|), our calculated value is 8 x 1CT 7 s _1 (cf. 
T able [2|, which is a factor of ~ 1.5 higher than that 
of Huebner et al. ( 1992 ) for solar minimum condi- 



tion at 1 AU. This difference in calculated values 
is mainly due to the higher (a factor of 1.24) value 
of solar flux at 1216 A in S2K model than that 
used by [Huebner et al~| ( |1992[ ) (cf. Figure [2]). 

No experimentally determined cross sections for 
the production of 0( 1 S) in photodissociation of 
H2O are available. The solar flux at H Lyman-a 
(cf. Fig. pi) is more than an order of magnitude 
larger than the flux at wavelengths below 1390 A, 
which is the threshold for the 0( 1 S) production in 
dissociation of H 2 0. To account for the production 
of 0( 1 S) in photodissociation of H2O, we assumed 
an yield of 0.5% at solar H Lyman-a (1216 A). 
However, to assess the impact of this assumption 
on the green and red line emissions we varied the 
yield between and 1%. The calculated photo- 
rate for the production of 0( 1 S) from H2O is 6.4 
x 10~ 8 s" 1 at 1 AU assuming 1% yield at 1216 A 
(cf. Table [l]). 

The primary dissociative product of H2O is 
OH. The important destruction mechanisms of 
OH molecule are pre-dissociation through fluores- 
cence process and direct photodissociation. The 
solar radiation shortward of 928 A can ionize OH 
molecule. The threshold wavelengths for the pro- 
duction of 0( 1 D) and 0( 1 S) in photodissociation 
of OH are 1940 and 1477 A, respectively. The dis- 
sociation channels of OH have been discussed by 



Budzien et al. ( 1994 1 and van Dishoeck & Dal- 
garno ( 1984[). We h ave used the photo-rates given 
by Huebner et al. (19921 for the production of 



0( 1 D) and 0( i S) from OH molecule whose values 
are 6.4 x 10~ 7 and 6.7 x 10 -8 s _1 , respectively. 
These rates are based on dissociation cross sec- 



tions of van Dishoeck & Dalgarno (1984), which 



are consistent with the red line observation made 



by wide-field spectrometer ( Morgenthaler et al. 



2007) 



3.1.2. C0 2 

The threshold wavelengths for dissociation of 
C0 2 molecule producing 0( 1 D) and 0( 1 S) are 



1671 A and 1286 A, re spectively. As noted by 
Huestis fc Slanger|(|2006), the 0( 1 D) yield in pho- 



todissociation of CO2 has never been measured 
because of the problem of rapid quenching of 
this metastable state. However, experiment by 



Kcdzierski et al. ( 1998 1 suggested that this disso- 



ciation channel can be studied in electron impact 
experiment using solid neon matrix as detector. 
Huebner et al.| ( 1992 ) estimated the cross section 
for 0( i D) production in photodissociative excita- 
tion of CO2 (see Figure Ej), and obtained photo- 
rate values of 9.24 x lO^'and 1.86 x 10~ 6 s" 1 for 
solar minimum and maximum conditions, respec- 
tively. Using S2K solar flux on 1996 March 30 our 
calculated rate for 0( 1 D) production in photodis- 
sociation of CO2 is 1.2 x 10~ 6 s _1 at 1 AU, which 



is higher than the solar minimum rate of Hueb 
ner et al~| ( |1992[ ) by a factor of 1.3. This variation 
is mainly due to the differences in the solar fluxes 
(cf. Figure[2]) in the wavelength region 950 to 1100 
A where the photodissociative cross section for the 
production of 0( 1 D) maximizes (cf. Figure B. 

|Lawrence| (19721 measured the 0( 1 S) yield in 
photodissociative excitation of CO2 from thresh- 
old (1286 A) to 800 A. The yield of 



et al. 



Lawrence 



( 1972 ) is dif ferent from that measured by Slanger 
(|1977[) in the 1060 to 1175 A region. How- 



ever, the yield from both experimental measure- 
ments closely matches in the 1110-1140 A wave- 
length region, where the yield is unity. In the ex- 



periment of Slanger et al. (1977), a dip in quan- 
tum yield is observed at 1089 A. Huestis et al. 



(2010) reviewed the experimental results and sug- 



gested the yield for 0( 1 S) in photodissociation 
of CO2. We calculated the cross section for the 
O^S) production in photodissociative excitation 
of CO2 (see Fi gure [4| by mul tiplying the yield rec- 
ommended by Huestis et al. (2010) with total ab- 
sorption cross section of C0 2 (see Figure[3|. Using 
this cross section and S2K solar flux, the rate for 
0( 1 S) production is 7.2 x 10~ 7 s" 1 at 1 AU. 

3.1.3. CO 

The threshold wavelength for the dissociation of 
CO molecule into neutral products in the ground 
state is 1117.8 A and in the metastable 0( 1 D) and 
C^D) is 863.4 A. Among the O-bearing species 
discussed in this paper, CO has the highest dis- 
sociation energy of 11.1 eV, while its ionization 



potential is 14 eV. Huebner et al. ( 1992 1 calcu 



lated cross sections for the photodissociative ex- 
citation of CO producing 0( 1 D) using branch- 
ing ratios from |McElroy fc McConnen] ( |1971| ) (cf. 
Fig. |4J) . Rates for the production of 0( i D) from 
CO molecule calculated by Huebner et al.| ( 1992 ) 
are 3.47 x 1(T 8 and 7.87 x 10 _b s _i for solar 
minimum and maximum conditions, respectively. 



Kedzierski et al. ( 1998 ) measured the cross sec- 



Using the cross section of Huebner et al. (1992) 



and S2K model solar flux, our calculated rate for 
the 0( 1 D) production from CO is 5.1 x 1CT 8 s' 1 
at 1 AU, which is 1.5 times higher than the so- 
lar minimum rate of Huebner et al. (1992). This 



difference in the calculated value is due to varia- 
tion in the solar fluxes used in the two studies in 
wavelength region 600 to 800 A (cf . Figure [2l . 

We did not find any reports on the cross section 
for the production of O^S) in photodissociation of 
the CO molecule. According to Huebner & Car- 
penter (1979) the rate for this reaction can not 
be more than 4 x 10~ 8 s _1 . We have used this 
value in our model calculations. This process can 
be an important source of 0( 1 S) since the comet 
Hyakutake has a higher CO abundance (~20%). 
Using this photorate and CO abundance, we will 
show that this reaction alone can contribute up to 
a maximum of 30% to the total 0( 1 S) production. 

3.2. Electron impact dissociation 

In our literature survey we could not find any 
reported cross section for the production of 0( 1 D) 
due to electron impact dissociation of H2O. |Jack-| 
man et al. ( 1977 ) have assembled the experimental 



and theoretical cross sections for electron impact 
on important atmospheric gases in a workable an- 
alytical form. The cross sections for electron im- 
pact on atomic oxygen given by |Jackman et al. 
( 1977 ) have been used to estimate emissions which 



leave the O atom in the metastable ( 1 D) state. 
The obtained ratios of 85% in ground and 15% 
in metastable state are used for the atomic states 
of C and O produced in electron impact disso- 
ciation of H2O, CO2, and CO. It may be noted 
that the ground state to metastable state produc- 
tion ratio of 89:11 is observed for atomic carbon 
and atomic oxygen produced from photodissocia- 
tion of CO ( Singh et al.|199l" ). However, as shown 
later, the contributions of these electron impact 
processes to the total production of 0( 1 D) are very 
small (<5%). 



tion for electron impact dissociative excitation of 
H 2 producing 0( 1 S), with overall uncertainty of 
30%. |LeClair fc McConkey| ( |1994[ ) measured cross 
section for the production of 0( 1 S) in dissociation 
of CO2 by electron impact; they claimed an uncer- 
tainty of 12% in their experimental cross section 
measurements. The cross section for fragmenta- 
tion of CO into metastable Of^S) atom by elec- 



tron impact is measured by LeClair et al. ( 1994 | 



These electron impact cross sections are also rec- 



ommended by McConkey et al. (2008), and are 



used in our model for calculating the production 
rate of 0( 1 S) from H 2 0, C0 2 , and CO. 

Since the *D and *S are metastable states, the 
direct excitation of atomic oxygen by solar ra- 
diation is not an effective excitation mechanism. 
However the electron impact excitation of atomic 
oxygen can populate these excited metastable 
states, which is a major source of airglow emis- 
sions in the upper atmospheres of Venus, Earth, 
and Mars. We calculated the excitation rates 
for these processes using electron impact cross 
sections from Jackman et al. (19771. In calcu- 



lating the photoelectron impact ionization rates 
of metastable oxygen states, we calculated the 
cross sections by changing the threshold energy 
parameter for ionization of neutral atomic oxy- 
gen in the analytical expression given by Jackman 
et al. ( 1977 ) . The above mentioned electron im- 



pact cross sections for the production of 0( 1 S) 
from H 2 0, CO2, CO, and O, used in the current 
model, are presented in Figure [5] along with the 
calculated photoelectron flux energy spectrum at 
cometocentric distance of 1000 km. 

3.3. Dissociative recombination 

The total dissociative recombination rate for 



H 2 0+ reported by |Rosen et al] ( |2000[ ) is 4.3 x 
10~ 7 cur 3 s- 1 at 300 K. The channels of disso- 
ciative recombination have also been studied by 
this group. It was found that the dissociation pro- 
cess is dominated by three-body breakup (H + H 
+ O) that occurs with a branching ratio of 0.71, 
while the fraction of two-body breakup (O + H 2 ) 
is 0.09, and the branching ratio for the formation 
of OH + H is 0.2. The maximum kinetic energy 
of the dissociative products forming atomic oxy- 
gen produced in ground state are 3.1 eV and 7.6 
eV for the three-and two-body dissociation, re- 



spectively. Since the excitation energy required 
for the formation of metastable 0( 1 S) is 4.19 eV, 
the three-body dissociation can not produce oxy- 
gen atoms in X S state. However, the 0( 1 D) atom 
can be produced in both, the three-body and the 
two-body, breakup dissociation processes. To in- 
corporate the contribution of ti20 + dissociative 
recombination in the production of O^D) and 
0( 1 S), we assumed that 50% of branching frac- 
tion of the total recombination in three-body and 
two-body breakups lead to the formation of 0( 1 D) 
and 0( 1 S) atoms, respectively. For dissociative re- 
combination of CO 2 , CO + and OH + ions we as- 
sumed that the recombination rates are same for 
the production of both 0( 1 D) and 0( 1 S). We will 
show that these assumptions affects the calculated 
0( 1 S) and 0( 1 D) densities only at larger (> 10 4 
km) cometocentric distances, but not in the inner 
coma. Tables [T] and [2] list the rates, along with 
the source reference, for these recombination reac- 
tions. 

4. Results and discussion 

4.1. Production and loss of 0( 1 S) atom 

The calculated 0( X S) production rate profiles 
for different processes in comet C/1996 B2 Hyaku- 
take are presented in Figure [6} These calculations 
are made under the assumption of 0.5% yield of 
0( 1 S) from H 2 at 1216 A solar H Lyman-a line 
and 1% CO2 relative abundance. The major pro- 
duction source of Of^S) is the photodissociative 
excitation of H2O throughout the cometary coma. 
However, very close to the nucleus, the photodis- 
sociative excitation of CO2 is also an equally im- 
portant process for the 0( X S) production. Above 
100 km, the photodissociative excitation of CO2 
and CO makes an equal contribution in the pro- 
duction of 0( 1 S). Since the cross section for elec- 
tron impact dissociative excitation of H 2 0, CO2, 
and CO are small (see Figure [5]), the contribu- 
tions from electron impact dissociation to 0( 1 S) 
production are smaller by an order of magnitude 
or more than that due to photodissociative excita- 
tion. At larger cometocentric distances (>2 x 10 3 
km), the dissociative recombination of H 2 + ion 
is a significant production mechanism for 0( S), 
whose contribution is higher than those from pho- 
todissociative excitation of CO2 and CO. The dis- 
sociative recombination of other ions do not make 



any significant contribution to the production of 
0( 1 S). 

In the inner coma, the calculated production 
rates of 0( 1 S) via photodissociative excitation is 
CO2 at various wavelengths are presented in Fig- 
ure Yn The major production of 0( 1 S) occurs 
in the wavelength region 955-1165 A where the 
average cross section is ^2 x 10~ 17 cm~ 2 (cf. 
Fig. H| and the average solar flux is ~1 x 10 9 
photons cm~ 2 s _1 (cf. Fig. pj. The calculated 
loss rate profiles of 0( 1 S) for major processes are 
presented in Figure [8] Close to the nucleus (<50 
km), quenching by H 2 is the main loss mecha- 
nism for metastable 0( 1 S). Above 100 km, the ra- 
diative decay of 0( 1 S) becomes the dominant loss 
process. The contributions from other loss pro- 
cesses are orders of magnitude smaller and hence 
are not shown in Figure [HJ 

4.2. Production and loss of 0( 1 D) atom 

The production rates as a function of cometo- 
centric distance for various excitation mechanisms 
of the 0( 1 D) are shown in Figure M The ma- 
jor source of O^D) production in the inner coma 
is photodissociation of H2O. The wavelength de- 
pendent production rates of O^D) from H2O are 
presented in Figure 10 The 0( X D) production in 



photodissociation of H 2 is governed by solar ra- 
diation at H Lyman-a (1216 A) wavelength. How- 
ever, very close to the nucleus, the production of 
0( 1 D) is largely due to photons in the wavelength 
region 1165-1375 A. Since the average absorption 
cross section of H2O decreases in this wavelength 
region by an order of magnitude, the optical depth 
at wavelengths greater than 1165 A is quite small 
(see Figure [3]) . Hence, these photons are able to 
travel deeper into the coma unattenuated, thereby 
reaching close to the nucleus where they dissoci- 
ate H 2 producing 0( 1 D). Thus, at the surface of 
cometary nucleus the production of 0( 1 D) is con- 
trolled by the solar radiation in this wavelength 
band. In high production rate comets, the pro- 
duction of 0( 1 D) near nucleus would be governed 
by solar photons in this wavelength region. The 
production of O^D) from H 2 by solar photons 
from other wavelength regions is smaller by more 
than an order of magnitude. 

After photodissociative excitation of H2O, the 
next significant O^D) production process at ra- 
dial distances below 50 km is the photodissociative 



excitation of CO2. Above 50 km to about 1000 
km, the radiative decay of 0( 1 S), and at radial 
distances above 1000 km the dissociative recom- 
bination of H2O" 1 ", are the next potential sources 
of the Of^D) (see Figure^]). The calculated wave- 
length dependent production rates of 0( 1 D) for 
photodissociation of CO2 are shown in Figure [Tl] 
Solar radiation in the wavelength region 1165-955 
A dominates the O^D) production. Since the 
cross section for the production of 0( 1 D) due to 
photodissociation of CO2 is more than an order of 
magnitude higher in this wavelength region com- 
pared to cross section at other wavelengths (see 
Figure El) , the solar radiation in this wavelength 
band mainly controls the formation of 0( 1 D) from 
CO2. Other potential contributions are made by 
solar photons in the wavelength band 1585-1375 A 
at distances <50 km, and 955-745 A at radial dis- 
tances >100 km. Since the C0 2 absorption cross 
section around 1216 A is smaller by more than 
two orders of magnitude compared to its maxi- 
mum value, the solar radiation at H Ly-a is not 
an efficient source of O^D) atoms. 



ties are presented in Figure 13 The 0( 1 D) density 



Zipf ( 1969 ) measured the total rate coefficient 
for the quenching of 0( 1 S) by H 2 as 3 x 10~ 10 
cm 3 s _1 . The primary channel in quenching mech- 
anism is the production of two OH atoms. The 
production of 0( 1 D) is also a possible channel 
whose rate coefficient is not reported in the lit- 
erature. Hence, we assumed that 1% of total rate 
coefficient can lead to the formation of 0( X D) in 
this quenching mechanism. However, this assump- 
tion has no implications on the O^D) production 
since the total contribution due to 0( 1 S) is about 
three orders of magnitude smaller than the major 
production process of 0( 1 D). 

The calculated loss rate profiles of 0( 1 D) are 
" Below 1000 km, the 0( 1 D) 
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presented in Figure 
can be quenched by various cometary species. The 
quenching by H 2 is the major loss mechanism for 
0( 1 D) below 500 km. Above 2 x 10 3 km radiative 
decay is the dominant loss process for 0( 1 D). 

4.3. Calculation of green and red-doublet 
emission intensity 

Using the calculated production and loss rates 
due to various processes mentioned above, and 
assuming photochemical equilibrium, we com- 
puted the number density of 0( 1 S) and 0( 1 D) 
metastable atoms. The calculated number densi- 



profile shows a broad peak around 200-600 km. 
But, in the case of 0( 1 S), the density peaks at 
much lower radial distances of ^60 km. The num- 
ber densities of 0( X D) and 0( 1 S) are converted 
into emission rate profiles for the red-doublet and 
green line emissions, respectively, by multiplying 
with Einstein transition probabilities as 



V(6300+6364)(>") = -4(6300+6364) X [O 1 D(r)} 



= A 



(6300+6364) ZZ^7 



T,i =1 Pi(r) 



YLiLi{r) + A?D 



(7) 



and 



V, 



r (5577)(r) = A {5577) x[0 1 S(r)} 



(8) 



Where [0 1 S{r)] and [0 1 D(r)] are the calculated 
number density for the corresponding production 
rates Pi(r) and loss frequencies Li{r) for 0( 1 S) 
and 0( 1 D), respectively. A( l D) and A^S) are the 
total Einstein spontaneous emission coefficients 
for red-doublet and green line emissions. Using 
the emission rate profiles, the line of sight inten- 
sity of green and red-doublet emissions along the 
projected distance z is calculated as 



I(z) 



n 



V, 



(5577, 6300+6364) 



(s)ds (9) 



where s is the abscissa along the line of sight, 
and V(5577 j 6300+6364) ( s ) is the emission rate for 
the green or red-doublet emission. The maximum 
limit of integration R is taken as 10 5 km. The 
calculated brightness profiles of 5577 and 6300 
A emissions are presented in Figure 
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These 

brightness profiles are then averaged over the pro- 
jected area corresponding to the slit dimension 
1.2" x 8.2" centred on the nucleus of comet 
C/1996 B2 Hyakutake for the observation on 30 
March 1986 ( |Cochran|20"08| . The G/R ratio aver- 
aged over the slit is also calculated. 

4.4. Model results 



Morrison et al. (1997) observed the green and 



red-doublet emissions on comet C/1996 B2 Hyaku- 
take in the high resolution optical spectra obtained 
on 1996 March 23 and 27 and found the G/R ratio 



in the range 0.12-0.16. Cochran (2008) observed 
the 5577 and 6300 A line emissions on this comet 
on 1996 March 9 and 30, with the G/R ratio as 
0.09 for 9 March observation. We calculated the 
G/R ratio by varying the yield for 0( 1 S) produc- 
tion in photodissociation of H2O at 1216 A (hence- 
forth refer to as 0( 1 S) yield). Since CO2 is not ob- 
served in this comet, we assumed that a minimum 
1% CO2 is present in the coma. However, we also 
carried out calculations for 0%, 3% and 5% CO2 
abundances in the comet. We calculated the con- 
tributions of different production processes in the 
formation of 0( 1 S) and O^D) at three different 
projected distances of 10 2 , 10 3 , and 10 4 km from 
the nucleus for the above mentioned CO2 abun- 
dances and the 0( 1 S) yield values varying from 
0% to 1%. These calculations are presented in 
Table [3] The percentage contribution of major 
production processes in the projected field of view 
for the green and red-doublet emissions are also 
calculated. The G/R ratio is calculated after av- 
eraging the intensity over the projected area 165 
x 1129 km which corresponds to the dimension 
of slit used in the observation made by |Cochran] 
(|2008h on 1996 March 30. These calculated values 



are presented in Table [4] 

Taking 1% C0 2 abundance and 0% 0( 1 S) yield, 
the calculated percentage contributions of major 
production processes of 0( 1 S) and 0( 1 D) atoms 
are presented in Table [3j Around 60 to 90% of 
the 0( 1 D) is produced from photodissociation of 
H2O. Contributions of photodissociative excita- 
tion of CO2 and CO in the production of 0( 1 S) 
and 0( 1 D) are 15 to 40% and 1%, respectively. 
Around 10 4 km projected distance, the photodis- 
sociative excitation of OH (~20%) and the disso- 
ciative recombination of H2 + (~30%) are also 
significant production processes for 0( S) atoms. 
But, the contributions from these processes in 
0( 1 D) production is around 10% only. 

For C0 2 abundance of 1% and 0( 1 S) yield 
of 0.2%, the calculations presented in Table [3] 
show that the photodissociation of H2O contribute 
around 20 to 40% in the production of 0( 1 S) and 
60 to 90% in the production of 0( 1 D) atom. The 
next major source of 0( 1 S) production is the pho- 
todissociation of CO2 and CO with each contribut- 
ing ~10 to 25%. The relative contributions from 
photodissociation of parent species H 2 0, CO2, and 
CO to 0( 1 S) and O^D) production decreases with 



increase in projected distance from the nucleus. 
At 10 4 km projected distance, the photodissoci- 
ation of OH contribute 15% and 8% to the pro- 
duction of 0( 1 S) and 0( 1 D) atoms, respectively. 
Above 1000 km projected distance, the contribu- 
tion of H20 + dissociative recombination to 0( 1 S) 
production is around 20%. The production of 
0( 1 D) atom is mainly via photodissociation of 
H 2 0, but around 10 4 km the dissociative recom- 
bination of H20 + ion is also a significant produc- 
tion process contributing around 12%. At 10 4 km, 
dissociative recombination of OH + also contribute 
around 10% to the total 0( 1 D) production, which 
is not shown in Table |3l and this value is indepen- 
dent of 0( 1 S) yield or CO2 abundance. Radiative 
decay of 0( 1 S) is a minor (<5%) production pro- 
cess in the formation of 0( 1 D). 

We also calculated the relative contributions of 
different processes in the formation of green and 
red line emissions in the slit projected field of view, 
which are presented in Table|4] For the above case, 
the photodissociation of H2O contribute around 
35%, while the photodissociation of CO2 and CO 
contribute 23% and 22%, respectively, to the pro- 
duction of green line emission. The contribution of 
dissociative recombination of H2 + ions is around 
10%. The major production process of red lines 
is photodissociation of H2O (90%); the dissocia- 
tive recombination of H 2 + and radiative decay 
of 0( X S) atom are minor (<5%) production pro- 
cesses. With the 0( 1 S) yield of 0.2% and 1% C0 2 
abundance, the slit-averaged G/R ratio is found 
to be 0.11. 

When the 0( 1 S) yield is increased to 0.5% with 
1% CO2 abundance (see Table[3]), the contribution 
from photodissociative excitation of H2O to the 
0( 1 S) production is increased, with value varying 
from 35 to 60%, while the contribution to 0( : D) 
production is not changed. In this case, the con- 
tribution from photodissociation of CO2 and CO 
to the O^S) production is reduced (values be- 
tween 10 to 15%). The contributions from other 
processes are not changed significantly. Table [4] 
shows that in this case around 60% of green line 
in the slit projected field of view is produced via 
photodissociation of H2O, while the contributions 
from photodissociation of CO2 and CO are around 
15% each. The main (90%) production of red- 
doublet emission is through photodissociation of 
H2O. The slit-averaged G/R ratio is 0.17. 
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On further increasing the 0( 1 S) yield to 1% 
with CO2 abundances of 1%, the contribution of 
photodissociation of H 2 to 0( 1 S) atom produc- 
tion is further increased (values between 50 to 
75%) while the contribution from photodissocia- 
tion of CO2 and CO is decreased to around 10% 
each (cf. Table [3]). The contributions from other 
processes are not affected compared to the previ- 
ous case. As seen from Table |4j in this case the 
contribution of photodissociation of H2O to green 
line is around 75% in the slit projected field of 
view, while contributions from photodissociation 
of CO2 and CO are decreased to 10% each. The 
calculated G/R ratio is 0.27 (Table |I]). 

We also evaluated the effect of CO2 on the red- 
doublet and green line emissions by varying its 
abundance to 0%, 3% and 5%. The calculated per- 
centage contribution of major processes along the 
projected distances and in the slit projected field 
of view are presented in Tables [3] and |4j respec- 
tively. In the absence of CO2, the contributions 
from H 2 0, H 2 0+ and CO in 0( 1 S) production 
are increased by ~10% (cf. Tables [3] and [1} . Tak- 
ing 0% 0( 1 S) yield and by increasing C0 2 relative 
abundance from 1 to 3%, the percentage contribu- 
tions for 0( 1 S) from photodissociative excitation 
of CO2 (CO) is increased (decreased) by 50%. The 
contribution from H2O to 0( : D) production is not 
changed. 

The calculations presented in Tables [3] and ^de- 
pict that the contributions of various processes are 
significant in the production 0( X S) atom, whereas 
photodissociative excitation of H2O is the main 
production process for 0( 2 D) atom. Since comet 
C/1996 B2 Hyakutake is rich in CO (abundance 
^22%) compared to other comets, the contribu- 
tion from CO photodissociation to 0( X S) produc- 
tion is significant (10-25%). In the case of a comet 
having CO abundance less than 20%, the ma- 
jor production source of metastable 0( X S) atom 
would be photodissociation of H2O and CO2. 

4.5. Comparison with observations 

In 1996 March, the green and red-doublet emis- 
sions were observed in comet C/1996 B2 Hyaku- 
take from two ground-based observatories ( Morri- 
son et al.fl997||Cochran||2008[ ). Each observatory 
determined the G/R ratio using different slit size. 
Using a circular slit, the projected radial distance 



tion on March 23 and March 27 varied from 640 to 



653 km, while for Cochran (2008) observation, us- 



ing a rectangular slit, the projected area was 480 x 
3720 km on March 9 and 165 x 1129 km on March 
30. The clear detection of both green and red- 
doublet emissions and determination of the G/R 
ratio could be done for March 9 and March 23 



observations only (Cochran 2008 Morrison et al 



1997). The observed G/R ratio was 0.09 and 0.12 



to 0.16 for the observation on March 9 and March 
23, respectively. 

Making a very high resolution (R = 200,000) 
observation of comet C/1996 B2 Hyakutake on 
1996 March 30, |Cochran| ( |2008[ ) obtained radial 
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profiles of 5577 and 6300 A lines. In Figure 
we have compared the model calculated intensity 
profiles of 6300 and 5577 A lines at different pro- 
jected cometocentric distances with the observa- 



tion of |Cochran| ( |2008[ ). The calculated G/R ra- 
tio along projected distance is shown in Figure [15) 
The 6300 A emission shows a flat profile upto ^500 
km, whereas the 5577 A green line starts falling off 
beyond 100 km. This is because of the quenching 
of 0( 1 S) and O^D) by H 2 in the inner most 
coma (cf. Figures [8] and 12), thereby making both 
the production and loss mechanisms being con- 
trolled by H2O. Above these distances, the emis- 
sions are mainly controlled by the radiative decay 
of *S and : D states of oxygen atoms. 

Similar to the calculations presented in Tables[3] 
and[4l in Figures 14 and 15 we present the red and 
green line intensity profiles and the G/R ratios, 
respectively, for different contributions of 0( 1 S) 
yield and C0 2 abundances. Since photodissocia- 
tive excitation of H 2 is the main production pro- 
cess for 0( X D) atom, the red line intensity is al- 
most independent of the variation in O^S) yield 
and CO2 abundance. In the case of 0% CO2 abun- 
dance, the best fit to the observed green line pro- 
file is obtained when the 0( 1 S) yield is ~0.5% (± 
0.1%), where the G/R ratio varied from 0.06 to 
0.26 (cf. Figure [l5| and the slit-averaged G/R ra- 
tio for March 30 observation is 0.15 (cf. Table|4|. 
The shape of green line profile cannot be explained 
with 1% or 0% 0( 1 S) yield, while the case for 0.2% 
0( X S) yield can be considered as somewhat consis- 
tent with the observation. For this case, the value 



over the comet for Morrison et al. (1997) observa- 



of G/R ratio shown in Figure 15 is found to vary 
over a large range of 0.54 to 0.02. 

When we consider 1% CO2 in the comet, the 
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best-fit green profile is obtained when the 0( 1 S) 
yield is ~0.2%. The case for 0.5% 0( X S) yield also 
provides the green line prohle consistent with the 
observation. In both these cases the G/R ratio 
varies between 0.32 and 0.04 over the cometocen- 
tric projected distances of 10 to 10 4 km. The cal- 
culated 5577 A profiles for 0( 1 S) yield of 0% and 
1% are inconsistent with the observed profile. 

In Figure [14] we also show a calculated profile 
for a case when the CO2 abundance is 3% while 
the 0( 1 S) yield is 0% (i.e., no 0( 1 S) is produced in 
photodissociation of H2O). The calculated 5577 A 
green line profile shows a good fit to the observed 
profile: suggesting that even a small abundance 
of CO2 is enough to produce the required 0( 1 S). 
This is because the CO2 is about an order of mag- 
nitude more efficient in producing 0( S) atom 
than H2O in the photodissociation process (see 
Table [l]). However, since 0( 1 S) would definitely 
be produced in the photodissociation of H2O, and 
that the CO2 would surely be present in comet 
(though in smaller abundance), the most consis- 
tent value for the 0( 1 S) yield would be around 
0.5%. Assuming 5% C0 2 and 0.5% 0( 1 S) yield, 
the calculated green line emission profile is incon- 
sistent with the observation (cf. Figure 15 1. In 



this case, the calculated G/R ratio shown in Fig- 
ure [15] is found to vary between 0.24 and 0.05. 

From the above calculations it is clear that the 
slit projected area on to the comet also plays an 
important role in deciding the G/R ratio. This 
point can be better understood from Table[5]where 
the G/R ratio is presented for a projected square 
slit on the comet at different geocentric distances. 
It is clear from this table that for a given physical 
condition of a comet and at a given heliocentric 
distance, the observed G/R ratio for a given slit 
size can vary according to the geocentric distance 
of the comet. For example, for a 0( 1 S) yield of 
0.2% (0.5%) and C0 2 abundance of 1%, the G/R 
ratio can be 0.17 (0.26) if the comet is very close to 
the Earth (0.1 AU), whereas the G/R ratio can be 
0.07 (0.1), 0.06 (0.08), or 0.06 (0.07), if the comet, 
at the time of observation, is at a larger distance 
of 0.5, 1, and 2 AU from the Earth, respectively. 
Further, a G/R ratio of ~0.1 can be obtained even 
for the 0( 1 S) yield of 0%. This suggests that the 
value of 0.1 for the G/R ratio is in no way a defini- 
tive benchmark value to conclude that H2O is the 
parent of atomic oxygen atom in the comet, since 



smaller (~5% relative to H2O) amounts of CO2 
and CO itself can produce enough O^S) compared 
to that from H 2 0. This table also shows that for 
observations made around a geocentric distance of 
1 AU, the G/R ratio would be generally closer to 
0.1. The G/R ratio observed in different comets 



ranges from 0.02 to 0.3 (e.g., Cochran 2008 Capria 
erat|[20T0l ). 



Thus, we can conclude that the G/R ratio not 
only depends on the production and loss mech- 
anisms of 0( X S) atom, but also depends on the 
nucleocentric slit projected area over the comet. 
Moreover, the CO2 plays an important role in the 
production of O^S), and thus the green line emis- 
sion, in comets. With the present model calcula- 
tions and based on the literature survey of dis- 
sociation channels of H2O, we suggest that the 
0( X S) yield from photodissociation of H2O can- 
not be more than 1% of the total absorption cross 
section of H2O at solar Ly-a radiation. The best 
fit value of 0( X S) yield derived from Figure[l4]for a 
smaller (1%) CO2 abundance in the comet C/1996 
B2 Hyakutake is 0.4 (±0.1)%. As per the Tables[l] 
and^l this means that the ratio of rates of 0( 1 S) 
to O^D) production in the H 2 photodissocia- 
tion should be 0.03 (±0.01), which is much smaller 
than the value of 0.1 generally used in literature 
based on Festou &; Feldman] ( 1981 1 . Further, if 
the source of red and green lines is CO2 (CO), the 
ratio of photorates for 0( 1 S) to 0( 2 D) would be 
around 0.6 (0.8) (see Tables [l]and[2]). 

To verify whether the 0( 1 S) yield of 0.5% (for 
the CO2 abundance of 1%) derived from Figure [T4J 
based on the comparison between model and ob- 
served red and green line radial profiles in comet 
Hyakutake on 1996 March 30, is consistent with 
the G/R ratio observed on other days on this 
comet, we present in Table [6] the G/R ratio calcu- 
lated for observations made on 1996 March 9, 23, 
27, and 30, along with the observed value of G/R 



ratio from Morrison et al. (1997) and Cochran 



(2008). These calculations are made by taking the 



solar flux on the day of observation using Tobiska 



(2004) SOLAR2000 model and scaled according 



to the heliocentric distance of the comet on that 
date. The CO abundance is 22%, same as in all 
the calculations presented in the paper. 

The calculated G/R ratio on March 9, when 
geocentric distance was 0.55 AU and H 2 pro- 
duction rate 5 x 10 28 s -1 , is 0.09 (see Table [6} 
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which is same as the observed ratio obtained by 
Cochranl d2008b. On March 23 and 27 the comet 



is closer to both Sun and Earth (geocentric dis- 
tance ~0.1 AU) and its H 2 production rate was 
4 times higher than the value on March 9. The 
calculated G/R ratio on March 23 is 0.12, which 
is in agreement with the observed ratio obtained 



by Morrison et al. (19971 



5. Conclusions 

The Green and red-doublet atomic oxygen 
emissions are observed in comet C/1996 B2 
Hyakutake in 1996 March when it was passing 
quite close to the Earth (A = 0.1 to 0.55 AU). A 
coupled chemistry-emission model has been devel- 
oped to study the production of green (5577 A) 
and red-doublet (6300 and 6364 A) emissions in 
comets. This model has been applied to comet 
Hyakutake and the results are compared with the 
observed radial profiles of 5577 and 6300 A line 
emissions and the green to red-doublet intensity 
ratio. The important results from the present 
model calculations can be summarized as follow- 
ing. It may be noted that some of these results 
enumerated below may vary for other comets hav- 
ing different gas production rate or heliocentric 
distance. 

1. The photodissociation of H2O is the domi- 
nant production process for the formation of 
0( 1 D) throughout the inner cometary coma. 
The solar H Ly-a (1216 A) flux mainly gov- 
erns the production of 0( 1 D) in the pho- 
todissociative excitation of H 2 0, but near 
the nucleus solar radiation in the wavelength 
band 1375-1165 A can control the formation 
of 0( 1 D) from H 2 0. 

2. Other than the photodissociation of H 2 
molecule, above cometocentric distance of 
100 km the radiative decay of 0( 1 S) to 
0( 1 D) (via 5577 A line emission), while 
above 1000 km the dissociative recombina- 
tion of H 2 + ions, are also significant source 
mechanisms for the formation of 0( 1 D) and 
0( 1 S) atoms. 

3. The collisional quenching of 0( 1 D) atoms by 
H 2 is significant up to radial distance of 
^1000 km; above this distance the radiative 
decay is the main loss mechanism of 0( X D) 



atoms. The collisional quenching of 0( X D) 
by other neutral species is an order of mag- 
nitude smaller. 

4. The photodissociation of H 2 is the major 
process for the production of 0( 1 S) atoms, 
but near the nucleus the photodissociation of 
C0 2 can be the dominant source. The solar 
H Ly-a (1216 A) flux controls the production 
of 0( X S) via photodissociative excitation of 
H 2 0. 

5. At cometocentric distances of <100 km, the 
main loss process for 0( 1 S) is quenching by 
H 2 molecule, while above 100 km the ra- 
diative decay is the dominant loss process. 

6. Since the photoabsorption cross section of 
C0 2 molecule is quite small at 1216 A, the 
contribution of C0 2 in the production of 
0( 1 S) and 0( X D) at the solar H Ly-a is in- 
significant. 

7. Because the C0 2 absorption cross section in 
the 1165-955 A wavelength range is higher 
by an order of magnitude compared to that 
at other wavelengths, the solar radiation in 
this wavelength region mainly controls the 
production of 0( 1 D) and 0( X S) in the pho- 
todissociative excitation of C0 2 . Moreover, 
the C0 2 absorption cross section in this 
band is also the largest compared to those 
of H 2 and CO. 

8. The cross section for the photodissociation 
of H 2 producing 0( 1 S) at the solar H Ly-a 
wavelength (with 1% 0( 1 S) yield) is smaller 
by more than two orders of magnitude than 
the cross section for the photodissociation 
of C0 2 producing 0( X S) in the wavelength 
region 1165-955 A . Though the solar flux 
at 1216 A is higher compared to that in the 
1165-955 A wavelength region by two orders 
of magnitude, the larger value of C0 2 cross 
section in this wavelength band enables C0 2 
to be an important source for the production 
of metastable 0( 1 S) atom. 

9. In the case of CO, the dissociation and ion- 
ization thresholds are close to each other. 
Hence, most of the solar radiation ionizes 
CO molecule rather than producing the 
0( 1 S) and 0( 1 D) atoms. 
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10. Though the CO abundance is relatively high 
(-22%) in comet C/1996 B2 Hyakutake, the 
contribution of CO photodissociation in the 
0( 1 D) production is small (~1%), while for 
the production of 0( 1 S) its contribution is 
10 to 25%. 

11. The photoelectron impact dissociative exci- 
tation of H2O, CO2, and CO makes only a 
minor contribution (<1%) in the formation 
of metastable 0( 1 S) and 0( 1 D) atoms in the 
inner coma. 

12. The 0( 1 S) density peaks at shorter radial 
distances than the 0( 1 D) density. The peak 
value of 0( 1 S) density is found around 60 
km from the nucleus, while for the 0(D) a 
broad peak around 200-600 km is observed. 

13. In a H20-dominated comet, the green line 
emission is mainly generated in the pho- 
todissociative excitation of H2O with con- 
tribution of 40 to 60% (varying according 
to the radial distance) to the total inten- 
sity, while the photodissociation of CO2 is 
the next potential source contributing 10 to 
40%. 

14. For the red line emission the major source 
is photodissociative excitation of H2O, with 
contribution varying from 60 to 90% depend- 
ing on the radial distance from the nucleus. 

15. The G/R ratio depends not only on the pro- 
duction and loss processes of the O^S) and 
0( 1 D) atoms, but also on the size of observ- 
ing slit and the geocentric distance of comet 
at the time of observation. 

16. For a fixed slit size, the calculated value of 
the G/R ratio is found to vary between 0.03 
and 0.5 depending on the geocentric distance 
of the comet. In the inner (<300 km) most 
part of the coma, the G/R ratio is always 
larger than 0.1, with values as high as 0.5. 
On the other hand, at cometocentric dis- 
tances larger than 1000 km the G/R ratio 
is always less than 0.1. 

17. The model calculated radial profiles of 6300 
and 5577 A lines are consistent with the ob- 
served profiles on comet C/1996 B2 Hyaku- 
take for 0( X S) yield of 0.4 (±0.1) and C0 2 
abundances of 1%. 



18. The model calculated G/R ratio on comet 
Hyakutake is in good agreement with the 
G/R ratio observed on two days in 1996 
March by two observatories using different 
slit sizes. 
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Table 1 
Reactions for the production and loss of 0( 1 S). 



Reaction 



Rate (cm 3 s 1 or s 



Reference 



H 2 
OH 4 
C0 2 
CO 4 
H 2 
OH 4 
C0 2 
CO 4 



► 0( 1 S) - 
O^S) + 

> o^s) - 

O^S) 4- 

©(iS) 4 
e ph -4> 0( 1 S) 
e ph -> 0( 1 S) 4 



h hi/ 
lw 

- hi/ 
hi/- 

h e p ; 



-H 2 

H 

- CO 
C 

4- others 

- others 
4- others 

others 



O + e ph -4> 0( 1 S) 



H 2 0+ 
OH+ - 

CO+ - 
CO+ - 

0( 1 S) 
0( 1 S) 
0( 1 S) 

oes) 

0( 1 S) 
0( 1 S) 
0( 1 S) 
0( 1 S) 
0( 1 S) 



2e 



0( 1 S) 
0( X S) 
0( 1 S) 
O^S) 
0+4- 

► 0+ - 
0( 3 P) 
OCD) 
^ 2 OH 

-> 0( 1 D) 

-> 0( 3 P) + 

-> 0( 1 D) 

> CO + 
-> 0( 1 D) - 

ew, -> 0( 1 D) + c 
-> 0( 3 P) + 
0^2 0( 1 D) 



6th 

etfc 

etfc 
h hi/ 
h e pft 



H 2 



C0 2 



CO 



- others 
others 
others 
others 



hi/ 2g72 

h^5577 

+ H 2 
C0 2 

+ co 2 

I- CO 



6.4 x 10" 
6.7x 10- 



(300/T e ) a5 
x (300/T e 
x (300/T e 
x (300/T e 



x 

.0.5 
,0.5 
.0.46 



7.2 x 
4.0 x 
9.0 x 
2.2 x 
4.4 x 

2.2 x 
3.0 x 

4.3 x 

6.3 x 
2.9 x 

5.0 x 
1.9 x 
2.7 x 
0.075 
1.26 

3 x 10~ 10 
3 x 10~ 10 

3.1 x 10 
2.0 x 10 
3.21 x 10~ 12 exp(-1327/T) 

7.4 x 10~ 14 exp(-961/T) 
8.56 x 10~ 9 

1.56 x 10~ 9 (T e /300) 094 
2.0 x 10 



lO" 7 

io- 8 

io- 10 

io- 10 

io- 8 

io- 10 

io- 8 

io- 7 

io- 9 

io- 7 

io- 8 

io- 7 

io- 7 



0.045 



x 0.01 



" n exp(-1330/T) 
11 exp(-1327/T) 



-14 



This work 



Huebner et al. 
This work 



(1992) 



Huebner & Carpenter ( 1979 1 

This work 

This work 

This work 

This work 

This work 



Rosen et al 



Gueberman 

Mitchell (1990) 




Mitchell 



(1990) 



Huebner et al 



1'his work 



(1992) 



Wiese et al. 
Zipf (1969) 



Atkinson "fc Welge 



Wiese et al. 



2IpT (1969) 



(1996) 



(1996) 



Capetanakis et al. 



Capetanakis et al". 



(1972 



1993) 



1993) 



Capetanakis et al. (1993) 
Berrington fc JJurkepjigsi I 
Berrington fc Burkej ffM I 



Krauss & Neumann (1975 



a This rate is calculated assuming 1% yield for the production of 0( 1 S) at 1216 A wavelength. 

b 0.045 is the assumed branching ratio for the formation of 0( 1 S) via dissociative recombination of 
H 2 0+ ion. 

c 0.01 is the assumed yield for the formation of 0( 1 D) via quenching of H 2 

Note. — The photorates and photoelectron impact rates are at 1 AU on 1996 March 30; e p h = photo- 
electron, e-th = thermal electron, hv — solar photon, T e = electron temperature, T = neutral temperature. 
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Table 2 
Reactions for the production and loss of 0( 1 D). 



Reaction 



Rate (cm 3 s l or s 



Reference 



H 2 + hv -> 0( 1 D) + H 2 
OH + \w -> 0( 1 D) + H 
C0 2 + hi/ -> 0( 1 D) + CO 
CO + hi/ -> 0( X D) + C 
O^S) -> O^D) + hu 557nm 
H 2 + 6p fc -> 0( 1 D) + H 2 + c 
OH + c ph -> 0( X D) + H + e 
C0 2 + e ph -> OfD) + CO + c 
C( X D) + c 



CO 
O 

H 2 0+ 
OH+ 

CO+ 
CO+ 

oes) 
o^s) 
oes) 
o^s) 

0( 1 D) 
0( 1 D) 
0( 1 D) 

o^d) 



t- e ph -»• 0( 1 D) 
e p/t -> 0( 1 D) 
+- e*fc - 
- e*/, -» 
eth -> 
e t h ->■ 
f e*h - 
f H 2 

+ co 2 
f co- 

+ hi/ - 



Qph 



OOD) + e«i 
0( 1 D) + H 2 



0( 1 D) + CO 
0( 1 D) + CO 



- 0( X D) + H 2 
0( 1 D) + H 
O^D) + CO 
O^D) + C( 1 D) 
0( 1 D) +e 
-> 0( X D) + H 2 
->■ 0( 1 D) + C0 2 

-> 0( 1 D) + CO 

^ 0+ + e 

0( 3 P)+ hj/eaoo 

0( 3 P)+ hi/ 6364 
-> 0+ + 2e 
-> 0( 3 P) + c 

> OH + OH 
-¥ 0( 3 P) + H 2 
->■ H 2 + 2 

> O + C0 2 
-> CO + 2 
+ CO 

^co 2 



8.0 x 10" 


-7 


6.4 x 10" 


-7 


1.2 x 10" 


-6 


5.1 x 10" 


-8 



8.1 x 
2.1 x 



io- 10 
io- 10 



0.5 



9.0 x 10" 



2.2 x 
7.4 x 



10~ 12 

io- 11 



2.0 x 10~ 10 
5.5 x 10~ 10 
8.0 x 10~ n 



(T e /300) 

exp(-120/T) 
exp(-625/T) 



This work 



Huebner et al. 
This work 
This work 



(1992) 



1.26 

2.1 x 10~ 10 

7 x 10~ n 

8.5 x 10~ 9 

7 x 10~ n 

3.7 x IO" 7 

4.3 x IO" 7 (300/T e ) 05 x 0.35 a 

6.3 x 10- 9 x (300/T e ) 048 
2.9 x IO" 7 (300/T e ) 05 

5 x IO- 8 (300/T e ) 046 
1.5 x IO- 10 (T e /300) 094 
3 x IO" 10 x 0.01 b 
2.0 x IO" 11 exp(-1327/T) 

7.4 x IO- 14 exp(-961/T) 
1.82 x IO" 7 

6.44 x IO" 3 
2.15 x IO- 3 
1.75 x IO" 7 



(1996) 



Wiese et al. 
This work 
This work 
This work 
This work 
This work 
Rosen et all (12000 
Cuebcrman (199 



Mitchell ( 1990 
Mitchell| fTM? 



BerringtorTSTBurke ( 1981 \ 
Zipf (1955J 



Capet anakis et al 



Capctanakis et al. 



Huebner et all T992 



Storey & Zeippen 
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Fig. 1. — Energy level diagram of atomic oxygen 
showing different spectroscopic transitions related 
to X S and 1 D states. 




200 400 600 800 1000 1200 1400 1600 1800 2000 
Wavelength (A) 

Fig. 2.— Solar EUV-UV flux from SOLAR2000 



(S2K) model of |Tobiska| p)04| for the day 1996 
March 30 at heliocentric distance of 1 AU. For 
compa rison the solar flux used by |Huebner et al.| 
(|1992[) is also shown. At H Lyman a (1216 A) 



the solar flux given by S2K model is higher than 
that of|Huebneretal.|(|1992|) by a factor 1.24. Sig- 



nificant differences in the two solar fluxes can be 
noticed in the wavelength ranges 800 to 1200 A, 
while above 1400 A the solar flux of S2K model is 



smaller than that of Huebner et al. ( 1992 1 




ami 

200 400 600 800 1000 1200 1400 1600 1800 2000 
Wavelength (A) 



Fig. 3. — The total photoabsorption cross sections 
for H2O, CO2 and CO are shown in the top panel 
and total photoionization cross sections are shown 
in the bottom panel. The cross sections are taken 



from Huebner et al. (1992|. 
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hv + CO, -> 0(D) 




200 400 600 800 1000 1200 1400 1600 1800 
Wavelength (A) 

Fig. 4. — Photodissociative excitation cross sec- 
tions for the production of 0( 1 D) from H 2 0, CO2, 
and CO. These cross sections are taken from lHueb^l 



ner et al. (1992). -A- represents the cross section 



value for the production of 0( 1 S) from H 2 at 
1216 A assuming 1% yield. 
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Fig. 5. — Cross sections for the production of 
metastable 0( 1 S) in electron impact dissociative 
excitation of H 2 0, CO2, and CO. Cross sections 
for electron impact excitation of O to Of^D) and 
0( 1 S) are also plotted. Calculated steady state 
photoelectron flux at 1000 km cometocentric dis- 
tance using S2K model solar flux on 1996 March 
30 is also shown with scale on right side y-axis. 
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Fig. 6. — Calculated radial profiles for major pro- 
duction mechanisms of 0( 1 S) along with the total 
production profile, hv — solar photon, e p h — pho- 
toelectron, and eth — thermal electron. 
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Fig. 7. — Calculated radial profiles for the pho- 
to-dissociation of CO2 producing 0( 1 S) at different 
wavelength bands. 
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Fig. 9. — Calculated radial profiles for the major 
production mechanisms of 0( 1 D) along with total 
0( 1 D) production rate profile. hv = solar photon, 
e ph = photoelectron, and e t h = thermal electron. 
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Fig. 8. — Calculated radial profiles for the major 
loss mechanisms of 0( 1 S) atom. 
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Fig. 10. — Calculated radial profiles for the pho- 
todissociation of H 2 producing 0( 1 D) at differ- 
ent wavelength bands. 
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Fig. 11. — Calculated radial profiles for the pho- 
todissociation of CO2 producing 0( X D) at differ- 
ent wavelength bands. 
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Fig. 13. — Calculated number density profiles of 
0( 1 S), 0( 1 D), 0( 3 P), and OH, along with those 
of H 2 0, CO, and C0 2 . 
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Fig. 12. — Calculated radial profiles for major loss 
mechanisms of 0( 1 D) atom. 
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Fig. 14. — Calculated emission brightness pro- 
files along projected distances for 5577 A (green) 
and 6300 A (red) line emissions for different 
CO2 relative abundance [CO2] and yield [Yld] for 
0( 1 S) production in photodissociation of H2O. 
The green and red emission intensities at different 
projected distances observed on March 30 taken 
from Cochran (2008) are also shown (filled sym- 



bols with error bars) for comparison with the cal- 
culated values. 
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Fig. 15. — Calculated green to red-doublet in- 
tensity ratio along projected distances for differ- 
ent CO2 relative abundance [CO2] and yield [Yld] 
for 0( S) production in photodissociation of H2O. 
G/R — calculated green to red-doublet intensity 
ratio averaged over slit projected size 165 x 1130 
km for C/1996 B2 Hyakutake on 1996 March 30. 
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